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Abstract 

If phase coherence determines the superconducting transition temperature 
T c in the cuprate oxides, it is of great interest to understand the role that 
dynamics of the phase fluctuation plays in bringing about depletion of the 
superconducting condensate. We will show that a phase correlation function 
can be calculated that allows us to describe depletion of superconducting con- 
densate as a result of the quantum fluctuation. In two dimension, dynamic 
phase fluctuation or pair fluctuation gives rise to a condensate depletion lin- 
ear with temperature as T— ► O in superconductors with nodes at the fermi 
surface. 

Superconducting order parameter has an amplitude and a phase. In the ground state 
the phase has to be blocked in all space. It has been suggested that the superconducting 
transition temperature T c in the oxide materials is related to the loss of long range phase 
coherence. Beyond T c the system lives, in this model, unlike classical superconductors 
in an ill understood mixture of Cooper pairs that have a nonzero amplitude and a highly 
fluctuating phase. The dynamics of the phase fluctuation and depletion of superfluid density 
of the superconducting state is the object of this study. 

We write the superconducting order parameter as ip (t) .where 

i/>{t) = \i/>\exp(iO(x,t)) (1) 
where \if)\ is the amplitude assumed constant in space and 6(x,t) is the phase assumed 



fluctuating in space and time.. The phase dynamics is contained in the simple hydrodynamic 
phase fluctuation hamiltonian 

H = ±«. (W) 2 + JL ( d p s ) 2 + \dp s e^ (t) (2) 

Here K s is the superfluid phase stiffness and Xp is the density susceptibility related to the 
density of states of paired electrons, dp s , a local density fluctuation that has long range 
coulomb interaction with nonlocal potential (p (t) and e*is the pair charge. We note that 9 
and the pairdensity p s are conjugate quantities so that all the dynamics of is contained 
in the excess pair charge e*dp s , the subscript s meaning superconducting phase. We have 
the Josephson relation governing the dynamics of the phase to induced local potential (t) 
through 

89 e* , , 

V 2 = e {2dp s - p qp ) (3) 

Here p qp is local quasiparticle population, if there is any . Let us assume that the coulomb 
interaction given by the third term of the hamiltonian (2) is zero. The resulting hamiltonian 
can be written as 

where p D is the condensate mass density, ia4s the superfluid flow, proportional to V6 1 (x, t) 
and c Q is the bare fermi velocity. The hamiltonian gives the elementary phase excitation as 
the Boguylubov- Anderson mode with the well-known dispersion w q = cq where is c =■%, in 
2d . In order to obtain the depletion of superfluid density due to dynamic phase fluctuation 
we use the usual trick of adding to the phase hamiltonian (4) an amplitude modulation term 
that represents an infinitesimal coupling to the medium , a source term 

#' = -^E^+^) ( 5 ) 

X 
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Where 77 is the interaction between pairs and the medium which in the end will be allowed 
to go to zero . The fluctuation of the longitudinal part of the order paramter is defined 
through 



A'(x) =A-A 

A = v + V>*;A = ^ + ^ 

a '--\/!e«, 2 



(6) 



Here n is the condensate density (n Q ~ ^ where is the hole density in oxide materials) , 
Q is an atomic volume in 3d (in 2dimension n Q is the pair density per unit area while Q is 
the area of an unit cell on the surface ), 6 q is the fourier transform of the phase fluctuation 
9. The total fluctuation hamiltonian reads 



H = H n + H 



This can be transformed in the fourier space as 
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The excitation spectra of the phase mode is significantly modified, due to coupling rj and is 
given by (after restoring h) 



w q (77) = hcq ( 1 



77m n Q 



h 2 n q 2 V f2 

At small momenta we now have a gap in the excitation spectra given by 



(9) 



w (r/) = c 
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(10) 



It is relatively straightforward to derive the relation for the condensate density n (T) at low 
temperature due to phase fluctuation. By writing the grand canonical hamiltonian 



H = H-fj,N 



(11) 



we can show that 



(i/>) = y/n {T)Sl 
Here G is the Gibb's free energy and is given by 



or) 

I r\=o 



dG (rj) 



(12) 



T)=0 



kTVt f 

G (77, T)=G (77, 0) + -j—^ J qdq In [1 - exp -(3w q (rj)] 



(13) 



This allows us to write in two dimension 



n Q (T) =n Q - v/n °^ f q dqN° 

^ J q 



drj 



(14) 



T)=0 



Here N° is the bose factor for the unperturbed energy given by 



N° = 

" exp(^) 



(15) 



Using equation (9) and substituting it into equation (14) , we obtain integrating equation 
(14) as T — > 0, the condensate density in two dimension 

n (T) — n 1 mkT 1 mkT 



n n 



27r n h 2 



(16) 



irN h h 2 

In this expression m is the electronic mass and N h is the total hole density. Care has to be 
exercised in using the expression (16) as n , the pair density at T=0 may be considerably 
less than Nh/2 due to quantum fluctuation which we calculate in the next section. In three 
dimension the same calculation gives a T 2 diminution of the condensate fraction as 0. 

The quantum phase fluctuation will be studied with an lattice hamiltonian equivalent 
to that of expression (4) that we have used earlier [4] .We have pointed out that almost all 
high T c cuprates are characterised universally as domains of coherence length £ separated by 
about the same distance from neighboring domains and that phase coherence between the 
domains is assured by condensate pair transfer between them at a rate given by the Josephson 
coupling energy J. This gives the discrete or coarse-grained version of the hamiltonian of (4) 
as a Josephson lattice hamiltonian 



H e = -jJ2(cos(6 l -6 j )) + 
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(17) 



The energy increase for a single cooper pair charge dp s = ±2e is given by in this pure 
transfer process 

_ (dp.f _ 2e* 

e '-^r~i; (18) 

e c is the quantum unit of energy increase for a single cooper pair transfer. At T = 0, the 
lattice hamiltonian is equivalent to a d + 1 system and is given exactly by 



H d e +1 = -Y. c ™^- e i) ( 19 ) 

y ij 

where a quantum coupling constant for quantum fluctuation can be defined given by 



9 = (20) 

We have shown [4] that beyond a critical coupling constant g = g c , quantum fluctuation 
destroys the superconducting ground state and drives the system insulating. We will derive 
a similar characteristic from the phase correlation function below. 
We define the phase correlation function at a point x by 

$ x = {exp[iO(t)-iO(0)]) (21) 

Since phase and charge are conjugate variables ,the charge operator p induces transition 
from the ground state \ to an excited state I with a probability given by the matrix 
element \{^e \p\ \I/g)| 2 suc1i that we have the sum rule 

J2(E-E G )\(V E \p\V G )\ 2 =e c (22) 

E 

The phase correlator can be written as 

exp-,^^) = $, (23) 

We will now calculate the phase correlator $ using the lattice hamiltonian (16). Consid- 
ering 6 as an operator we write the commutator of 6 (0) and 9 (t) 



[9(Q),9(t)] = —anw t (24) 
Xp 

where w a is a characteristic frequency of the phase mode. We now use the operator 
identity 

A B [AM 



exp(A + B) = e A e B e 



to obtain 



exp id (0) exp — i9 (t) 

= exp i [9 (0) - 9 (t)) exp ( sin w t ) (25) 

V J 

we note that at small time 9 (t) ~ 9 (o) and the correlator behaves like exp (i^) as it 
should. 

This expression can be fed back into (^3|) to get after a little algebra 

/ (E — E )t\ n=+oo 

/exp— i — - — \ = $ x = P(nhw ) exp —inzu t (26) 

' ' n=— oo 

The fourier transform P(nfra7 ) gives the probability that n quanta of phase oscillation 
mode w Q will be excited in the medium. This is explicitly written in terms of the Bessel 
functions I n 

PC*,.) = exp - I„ ( ] —± m ) exp - (^) (27, 

Here exp — (6* 2 ) is like the Debye-waller factor and is given by the ratio of two charac- 
teristic enegies of the problem 

W = ife = ' < 28 ' 

where a quantum coupling constant has been defined which tells us how many quanta 
will be excited during a typical pair transfer so as to cause dephasing and depletion of 
the superconducting condensate. The expression for P(0) gives us the probability of 
condensate fraction at a given T. 
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This allows us to write for the spectral density at T=0 as 

( n=co 9 n \ 
A(u) = (exp -g) ^ {u-e c - nw Q ) (29) 

\n=0 U - ) 

.The spectral density obey the following sum rules 

> 

ujA (u) duj = e c (30) 

u m A {uj) du = (u m ) (31) 

Expression fl29|) is plotted in figure 1 for a variety of values of g. The rapid depletion of 
spectral weight A (u — o) i.e that of superfluid density at around g=g c is to be seen. The 
two quantum coupling constants of the expressions (28) and (20) are the same if we identify 
the characteristic energy t£7 with the average energy scale of the Josephson hamiltonian (19) 
as zu = \J t c .l\ w Q is just then a debye energy or a characteristic pair fluctuation frequency. 
We notice from the figures how the spectral weight goes from a poisson distribution to a 
gaussian as the coupling constant increases. The gaussian spectral distribution indicates that 
the ground state is coupled to a bath of continuum of excitations. We also note that the 
spread of the spectral weight as g increases does not mean we are in the excited state for 
we are still at T=0. If there were no interactions (g=0) we will get a single delta function 
at uj = e c . For g^ we are still in the ground state but the coupled system has n quanta of 
phase oscillation modes excited with a certain probability P(nhw Q ) .The effective superfluid 
density at T=0 can be written as 

■Sr ~ exp -g (32) 

We show in table 1 this fraction for a variety of optimum doped materials taken from 
||.In the table the effective number of holes per copper, n e ff has been measured by summing 
the oscillator strength N e //(o;) upto a frequency of the order of charge transfer band (~ 12000 
cm -1 ). As temperature is reduced, the increase in area at low frequencies is compensated by 
a decrease in area at high frequencies. However, below T c ,there is a loss of area, which must 



be coming mostly from the free carrier peak at to — > 0, as the authors conclude. Tne number 
of supeconducting carriers n s is obtained in two ways: from the London penetration depth 
Xl and from the Drude peak, in the normal state. The measured n s values are almost same 
as the number of Drude carriers nd(these later are not shown here). Thus only a quarter to 
a fifth of the total doping induced holes appear in the delta function !A value of g~ 1 can 
be taken indicating that even at optimum doping considerable quantum fluctuation is going 
on. It is interesting to compare this ratio of ^^with Ornstein et al's || optical data where 
a bulk plasma frequency from band structure value was estimated to be ~ 3.5e.v compared 
to 1.4 e.v of the superfluid condensate at uj=0 for YBa Cuo and which reflects faithfully the 
ratio observed in the table 1. 



Tablel 



material 




n e// 
CU 


n s 

CU 


n s % 

n eff 


La2cuc>4 


40 


0.15 


0.028 


19 


BiPbSrc«208 


79 


0.37 


0.085 


23 


BiSrCacii 2 Og (a) 


85 


0.44 


0.100 


23 


BiSrCacM20g (b) 


85 


0.48 


0.090 


19 


YBa 2 cu 3 07 (a) 


91 


0.44 


0.096 


22 


YBa 2 CM 3 07 (b) 


91 


0.59 


0.125 


21 


T\ 2 Ba 2 cacu20s 


110 


0.54 


0.115 


21 



The underlying physics of quantum phase fluctuation is quite similar to the boson shake 
up process in the strong coupling polaron problem |J 

We are now in a position to calculate the temperature coefficient of depletion of superfluid 
density using expression(16) and (33) and renormalising the temperature T by r = ^- so 
that we get the depletion coefficient in the dimensionless form as 

d {^) exp(,)^T e 
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Using for high T c supercoductors the following numbers T c = 100°k,m = 5 electron mass, 
Nh ~ 10 14 /cm 2 and g = 1, we obtain the coefficient as ~ — 0.56 in reasonable agreement 
with the values in the literature 0. The ratio on the right hand side is mkT c /n s , which 
is really a constant as shown by Uemura's plot. A linear T decrease of superconducting 
condensate due to static phase fluctuation has been suggested in the past || and that has 
been recently contested || .There is also the theory of low energy quasi-particle excitation 
HI at the d-wave gap nodes explaining the linear T-dependence but the numerical agreement 



seems to need considerable fermi-liquid correction || to fit the data. 

The overall treatment in this paper has assumed perfect screening of charge fluctuation 
dp s in the condensate so that plasma mode is not excited, we can speculate that in d-wave 
superconductor this may well be so due to possibility of low energy quasiparticle excitation 
as a result of phase excitation. The latter shifts the quasiparticle energy due to doppler 



effect [ll]of the local superfluid flow by an amount 



dE qp = hk~fv^ 

at the fermi surface i.e. at the nodal points. Here v s = — (h/2m) V0, is the local 
superfluid velocity. The quasiparticle population created due to this doppler shift is 
Pqp ~ f® Eqp f(E)N(E)dE where N(E), the quasiparticle density of states. The inte- 
gral picks up significant weight only because of d-wave nodes, where the product N(E) 
f(E) is 7^ 0. The condition of perfect screening (local charge neutrality) p qp « dp s may be 
easily satisfied so that long range coulomb interaction is suppressed. A background quasi- 
particle density of the order of 10% has been inferred from microwave measurements on 
high T c oxides ||12||.The local charge neutrality condition implies div J s +div J qp = , where 



J's are superconducting and quasiparticle currents respectively. For intense local phase gra- 
dient , the more stringent well-known criterion of total ( although local), destruction of 
condensate, due to quasiparticle backflow current applies with Jf lt = -J^ rf . 
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Figure 1: spectral density A (u>) as a function of phase quanta n. a) coupling constant 
9 = 1; b) g = 3; c) g = 8. 
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